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Abstract 

We present the first calculations to follow the evolution of all stable isotopes (and their 
abundant radioactive progenitors) in a finely zoned stellar model computed from the onset 
of central hydrogen burning through explosion as a Type II supernova. The calculations 
were performed for a 15 Mq Pop I star using the most recently available set of experimental 
and theoretical nuclear data, revised opacity tables, and taking into account mass loss due 
to stellar winds. We find the approximately solar production of proton-rich isotopes above 
a mass number of A = 120 due to the 7-process. We also find a weak s-process, which 
along with the 7-process and explosive helium and carbon burning, produces nearly solar 
abundances of almost all nuclei from A = 60 to 85. A few modifications of the abundances 
of heavy nuclei above mass 90 by the s-process are also noted and discussed. New weak rates 
lead to significant alteration of the properties of the presupernova core. 



1.1 Introduction 

Stars above ~ 10 M Q are responsible for producing most of the oxygen and heavier elements 
found in nature. Numerous studies of such stars and their detailed nucleosynthetic yields for 



various masses and metallicities, have been carried out previously, e.g., [32, 27]. However, our 
knowledge of both the input data and the physical processes affecting the evolution of these 
stars has improved dramatically in recent years. Updated opacity tables Jl^] have become 
available along with more accurate prescriptions for mass loss due to winds and new weak 
rates [|l(| that affect the evolutionary stages after central oxygen depletion. Perhaps most 
important for nucleosynthesis, new, accurate reaction rates for all the relevant strong and 
electromagnetic nuclear reactions above neon have been recently calculated by Rauscher and 



Thielemann [21|. Here we present the first results for a 15 M supernova evolved with the 
new physics. Additional masses and metallicities will be explored in future papers. These 
future papers will also include rotationally induced mixing processes ||. 

We also employ a nuclear reaction network of unprecedented size. The nuclear reaction 
network used by |32| (WW95), large in its day, was limited to 200 isotopes and extended 
only to to germanium (see also Chieffi and Limongi in these proceedings). Studies using 
nuclear reaction networks of over 5000 isotopes have been carried out for single zones or 
regions of stars, especially to obtain the r-process, e.g., [|], |5|, but "kilo-isotope" studies of 
nucleosynthesis in complete stellar models (typically of 1000 zones each) have been hitherto 



2 



lacking. We thus also present the first calculation to determine, self-consistently, the the 
complete synthesis of all stable isotopes in any model for a massive star. However, because 
its thermodynamic properties continue to be poorly determined (and for lack of space), we will 
ignore here the nucleosynthesis that occurs in the neutrino wind, which may be the principal 
site of the r-process |3jJ . 



1.2 Input Physics 

Our calculations were performed using a modified version of the stellar evolution code KE- 
PLER E32| with the following modifications: 



• updated neutrino loss rates ]1 

• improved opacity tables (OPAL95) g|, HI 

• mass loss due to stellar wind IT 

• updated weak rates fl 



updated strong and electromagnetic reaction rates |21] 



As in WW95, nucleosynthesis was followed by co-processing the stellar model throughout 
its evolution using an extended nuclear reaction network. From hydrogen ignition through 
central helium depletion a 617 isotope network was employed that included all elements up to 
polonium, adequate to follow the s-process. Just prior to central carbon ignition, we switched 
to a 1482 isotope network (also including astatine). This new network incorporated more 
neutron-rich isotopes to follow the high neutron fluxes in carbon (shell) burning. Five to 
ten isotopes were also added on the proton-rich side, in particular for the heavy elements, to 
follow the 7-process |3(| |2^, The nucleosynthesis during the supernova explosion itself 
was followed in each zone using a 2439 isotope network that included additional proton-rich 
isotopes to follow better the 7-process in the neon-oxygen core, and also many additional 
neutron-rich isotopes to follow the n-process expected during supernova shock front passage 
through the helium shell ||. 

We implemented a new library of experimental and theoretical reaction rates. In particu- 
lar, we used theoretical strong rates from [21] (using input from the FRDM |l7j| ), experimental 



neutron capture rates (30 keV Maxwellian average) along the line of stability by |TJ , and ex- 
perimental and theoretical rates for elements below neon as described in 111]]. Experimental 
(a, 7) rates were implemented for 70 Ge || and 144 Sm [25|. The derived a+ 70 Ge and a+ 144 Sm 



potentials were also utilized to recalculate the transfer reactions involving these potentials. 

Experimental fi~ , f3 + , and a-decay rates were taken from [^TJ, experimental (3~ rates from 
|l4|| and ||26|| , and theoretical f3~ and (3 + rates from ]l8|. As a special case, we implemented 
the temperature-dependent 180 Ta decay as described in 0. 



The supernova explosion was simulated, as in [32 1, by a piston that first moved inward for 
0.45 s down to a radius of 500 km and then moved outward to a radius of 10 000 km such that a 
total kinetic energy of the ejecta at infinity of 1.2- 10 51 erg resulted (for the 25 M stars we used 
a total kinetic energy of 1.5- 10 51 erg). The final mass cut outside the piston was determined 
by the mass that had settled on the piston at 2.5-10 4 s after core collapse. Note that the 
amount of fallback resulting from this prescription depends on both the initial location of the 
piston used and the energy of the explosion. In particular, the yields of 44 Ti and 56 Ni are 
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very sensitive to this "final mass cut" determined by the fall back. Multidimensional effects of 
the explosion are not considered here. The temperature of the \i and r neutrinos emanating 
from the proto-neutron star and causing the z^-process nucleosynthesis [34] were assumed to 
be 6MeV in contrast to WW95 who assumed 8MeV. However, we do not follow the ^-process 
for isotopes with Z or N larger than 40. 



1.3 Results and Discussion 
1.3.1 Stellar Structure 

Table Q gives the presupernova properties of our new models and, for comparison, those of 
WW95. The helium, carbon-oxygen, and neon-oxygen cores were defined as the location 
where hydrogen, helium, and carbon, respectively, first drop below a mass fraction of 1%, 
from the stellar surface going inward. The silicon core was defined by where silicon becomes 
more abundant than oxygen and the iron ( "Fe" ) core by where the sum of the mass fractions 
of 48 Ca and heavier nuclei first exceeds 50 %. The deleptonized core was defined as the region 
where the number of electrons per baryon, Y e , first drops below 0.49. 

The lower helium core masses in the new models (Table |]) are due to both mass loss and 
the use of the OPAL opacities. In the 25 M case a model using OPAL opacities but no mass 
loss resulted in a helium core of 8.69 M . As a result of the reduced helium core size our new 
models typically have lower helium-free and carbon-free cores. Due to the interaction of the 



Table 1: Properties of stellar models at the onset of core collapse (first section) and integrated 
stellar yields of some important radioactive nuclei (second section). 
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Figure 1: Production factors of iron group and lighter nuclei in a 15 M star of solar metal- 
licity. Shown here are the average integrated abundances in the ejecta (including mass loss 
by stellar winds) relative to solar Q (production factors). The dashed line indicates the pro- 
duction factor of 16 and the dotted lines span a band of ±0.3 dex. 



different phases of shells burning, the sizes of the "inner cores" do not always monotonically 
change with the size of the helium core. More details will be given elsewhere [22, 1. 

Some important changes in the new models are due to the revised weak rates | lq| . These 
rates become important during core silicon burning and thereafter. Typically, they lead to 
an increase of the central Y e at the onset of core collapse by 2 to 3 % (Table [j]) , and this 
difference tends to increase with increasing stellar mass [||. Perhaps more important for the 
explosion mechanism of core collapse supernovae is an increase of the density in the mass 
range of m = 1.5 M to 2M by 30 — 50% relative to the same models computed with the 
previous set of weak rates [32]. This may significantly affect the dynamics of the core collapse. 
For further details see 



1.3.2 Nucleosynthetic Production Factors 

In Figs. |] through ^] we show the production factors of all ejecta of the star after the explosion, 
including all the mass lost due to stellar winds, relative to solar Q abundances. We assume 
that all radioactive nuclei have decayed to their stable products. As a gauge we provide the 
production factor of 16 0, the dominant "metal" produced in massive stars (dashed line), and 
a band of acceptable agreement of ±0.3 dex relative to this values (dotted lines). 
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Figure 2: Production factors of trans-iron group nuclei (all isotopes are on scale) 
1.3.3 Light Elements and the Iron Group 

The species 2 H, 3 He, lithium, beryllium and boron were destroyed in the envelope of the star 
during central hydrogen burning. However, substantial 7 Li and n B were recreated by the 
^-process during the explosion [34], as was 19 F (Fig. [l]). 17 was significantly underproduced 



as a result of the revised reaction rates for 17 0(p, a) 14 N and 17 0(p,7) 18 F [|j 

The isotopes 18 through 38 Ar are in good agreement with solar abundance ratios. 40 Ar 
and 40 K are both significantly higher while other potassium isotopes are lower. This signature 
for the potassium was also found in other stellar models. The under-abundance of 44 Ca was 
caused by the low yield of 44 Ti (Table [fj) which beta-decays to calcium. The yield of this 
isotope strongly depends on the location of the final mass cut, i.e., the amount of fall back, 
and might also be affected by mixing processes during the supernova explosion. The same 
caveat also applies for other isotopes that mainly originate from regions close to the neutron 
star, like 56 Ni (see contribution of Kifonidis in these proceedings). 

1.3.4 The s-Process 

The nuclei above the iron group up to about A = 90 (Fig. ^) are produced as secondary 
isotopes by the s-process starting from iron. When considering galactic chemical evolution 
these yields are to be combined with those of metal-poor stars that contribute correspondingly 
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Figure 3: Post iron group nuclei 

less of these isotopes, therefore a production factor of about twice that of 16 is in good 
agreement with reproducing the solar abundance pattern. Note that the yields of these 
isotopes, by abundance, starting from iron decreases about exponentially. 64 Zn, which is 
underproduced as shown in Fig. |2[ may be a product of the neutrino wind from the proto- 
neutron star JlO[ ] . The possible contributions due to this process are not included in the 
results presented here. 

The overabundance of the neutron-rich nickel isotopes, 61 ' 62 ' 64 Ni, and other s-process 
products in the A = 60 — 90 mass range has been observed before |28|, |ll|] and is still not well 
understood. It is even greater in stars of 20 and 25 M©. Perhaps the problem will be alleviated 
by a more complete grid of supernovae of various metallicities and masses, perhaps the stellar 
structure will be altered by still uncertain physics (overshoot, 12 C(a,7) 16 0, rotation), or 
perhaps key reaction rates responsible for neutron production or absorption will change. For 
now, it remains problematic. 

Above A = 100 (Figs. || and |j) the s-process had only minor effects in this 15 M© star, 
though there were important "redistributions" of some of the heavy isotopes. Most of the 
s-process above mass 90 is believed to come from AGB stars. 

The "cutoff" towards lower values at a production factor of ~ 0.9 is due to the fact that 
most of the star does not become hot enough to affect the abundances of these nuclei, or is 
even lost in the wind. The supernova ejecta containing regions depleted by the s-process and 
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Figure 4: Heavy nuclei 

other processes are then averaged with the dominating contribution of unaffected matter. In 
the 15 M star presented here, this leads to the fact that 80% of all ejecta, including winds, 
did not experience the s-process. 



1.3.5 The 7-Process 

The production of the proton-rich nuclei results from photo-disintegration of heavy nuclei 
during implosive and explosive oxygen and neon burning (7-process [30, |2^, ^ 



|2Jj). Here we 

present the results of the first calculations that follow the 7-processes through the presuper- 
nova stages and the supernova explosion in the whole star. Fully self-consistent, the 7-process 
here operates in stellar regions that were exposed to previous episodes of s-processing. 

Above A = 123 to A = 150 and between A = 172 and A = 200 the proton-rich heavy 
isotopes are produced in solar abundance ratios within about a factor of two relative to 16 
(Figs. ^ and||). Below A = 123 and around A = 160 the production of the proton-rich isotopes 
is down by about a factor of three to four. The total production of the proton-rich isotopes 
increases for higher entropy in the oxygen shell, i.e., with increasing mass of the helium core, 
as we have seen in our 25 M star, but also depends on details of stellar structure and the 
composition of the star at the time of core collapse. Therefore the contribution from more 
massive stars may well fill in the gaps of the low production factors seen in the 15 M Q star. 

The isotope 180 Ta, the rarest stable nuclear species in the solar abundance pattern, shows 
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a remarkable overproduction (Table ||) in all of our models, despite our taking into account 
its destruction by de-excitation into the short-lived ground state through thermal excitation 
into an intermediate state P]. This may indicate that decay from other excited states could 
be important, which are not accounted for here. We cannot exclude, however, that our 
treatment of 180 Ta as a single species in the excited state only may cause, at least in part, 
the overproduction found here. 



1.3.6 The r- and n-Process 

The base of the helium shell is suspected to be a possible site for fast neutron capture processes 
as the supernova shock front passes these layers, especially in the less massive core collapse 
supernovae. Since current models of r-process sites have difficulties in reproducing the r- 
process peak around A = 130 when adjusted to fit the heavier nuclei (see the contribution 
of Truran in these proceedings), the base of the helium shell was considered as a possible 
environment for producing these isotopes. 

In our present models a distinct redistribution of nuclei around A = 123 was found at 
the base of the helium shell, but the resulting yields were too small to constitute a signifi- 
cant contribution that would be visible in Fig. [3|. We may speculate that less massive core 
collapse supernovae might have a stronger contribution though. More details on the present 



calculations will be given in [22] 



1.4 Conclusions 

We have presented the first calculation to follow the complete s-process through all phases 
of stellar evolution and the 7-process in the whole star through the presupernova stage and 
subsequent supernova explosion. Below, we summarize the important results for our 15 Mq 
star. Note, however, that though this mass is a numerically typical case of a Type II or lb 
supernova, the average nucleosynthetic yield of massive stars is the result of populations of 
different stars each of which has its own peculiar yields which must be combined to result in 
a solar-like abundance pattern. Some isotopes that are underproduced here may be strongly 
overproduced in other massive stars while isotopes overproduced here may be deficient in 
others. 

The proton-rich heavy isotopes above A = 123 can be well produced by the 7-process 
occurring during implosive and explosive oxygen and neon burning. The proton-rich isotopes 
around A = 160 and those between A = 100 and A = 123, however, are underproduced by a 
factor of 3 to 4 with respect to 16 0. The isotope 180 Ta shows a strong overproduction by the 
7-process. This may indicate that decay from excited states, of which we include only one, 
could be important. 

A strong secondary s-process contribution appears between iron and a mass number of 
A = 90. Above A = 100 the s-process in our 15 M star is very weak, but it becomes notably 
stronger in stars with more massive helium cores that perform helium burning at higher 
entropies. 

The expected r- or n-process contribution due to the supernova shock front running 
through the base of the helium shell does not show a significant contribution in any of our 
preliminary model stars, not even at A = 130. We observed some redistribution of isotopes 
at the base of helium shell around A = 123, but this did not show the characteristics of a 
typical r-process nor was it important compared to the total yield of the star. 
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The revisions of opacity tables and the introduction of mass loss generally leads to smaller 
helium core sizes which tend to also decrease the mass of the carbon-oxygen and the silicon 
core (Table ||). Note, however, that the absolute values of these core masses depend on the 
uncertainties, in particular, of the mixing processes in the stellar interior, such as semicon- 
vection, overshooting, and rotation. 



The revision of the weak rates [16], important after central oxygen burning, leads to a 
2 — 3 % higher electron fraction per nucleon, Y c , at the time of core collapse in the center of the 
star (Table |l|) and the "deleptonized core" tends to comprise less mass [||. More important 
for the core collapse supernova mechanism might be the 30 — 50% higher densities of the 
new models between the region of m = 1.5 — 2 M Q ||], which may result in a correspondingly 
higher ram-pressure of the infalling matter. 
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